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1. Introduction 

One of the most spectacular predictions of general relativity are black holes, which 
should form when a large mass is concentrated in a sufficiently small volume. The idea 
of a mass-concentration, which is so dense that even light would be trapped goes back 
to Laplace in the 18th century. Shortly after Einstein developed general relativity, Karl 
Schwarzschild discovered in 1916 a mathematical solution to the equations of the theory 
that describes such an object. It was only much later, with the work of people like 
Oppenheimer, Volkoff and Snyder in the 1930’s, that the scientific community began 
to think seriously about the possibility that such objects might actually exist in the 
Universe. It was shown that when a sufficiently massive star runs out of fuel, it is 
unable to support itself against its own gravitational attraction and it should collapse 
into a black hole. Only in the 1960’s and the 1970’s, in the so-called Golden Era of black 
hole research, new interesting phenomena like the Hawking radiation and superradiance 
were discovered but for a rigorous mathematical description of them we have to wait 
until the 1990’s and the beginning of the new century, when the rigorous analysis of the 
propagation and of the scattering properties of classical and quantum fields on black 
hole space-times was developed. 

Whenever we attempt to analyze the scattering properties of fields in the more 
realistic framework of the Kerr-Newman black hole geometry, we are faced with several 
difficulties, which are not present in the simplified picture of the Schwarzschild metric. 
First of all, the Kerr-Newman solution is only axisymmetric (cylindrical symmetry), 
since it possesses only two commuting Killing vector fields, namely the time coordinate 
vector held d t and the longitude coordinate vector held d This implies that there is 
no decomposition in spin-weighted spherical harmonics and that artihcial long-range 
terms at inhnity will be present in the held equations. Moreover, another difficulty 
is due to fact that the Kerr-Newman space-time is not stationary. In particular it is 
impossible to fold a Killing vector held which is time-like everywhere outside the black 
hole. In fact dt becomes space-like in the ergosphere, a toroidal region around the 
horizon. This implies that for held equations describing particles of integral spin (wave 
equation, Klein-Gordon, Maxwell) there exists no positive definite conserved energy. 
For held equations describing particles with half integral spin (Weyl, Dirac) we can 
fold a conserved L 2 norm with the usual interpretation of a conserved charge. Hence, 
the absence of stationarity in the Kerr-Newman metric is not really a difficulty for the 
scattering theory of classical Dirac helds. 

For the reasons mentioned above there are only few analytical studies of the propagation 
of helds outside a Kerr-Newman black hole. Time-dependent scattering for the Klein- 
Gordon equation in the Kerr framework has been developed in 2003 by Hafner T], while 
in 2004 Daude [2] proved the existence and asymptotic completeness of wave operators, 
classical at the event horizon and Dollard-modihed at inhnity, for classical massive Dirac 
particles in the Kerr-Newman geometry. 

Our present work is the hrst in a series of three papers devoted to develop a time- 
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dependent scattering theory for massive Dirac particles outside a non-extreme Kerr- 
Newman background. Our main technical tool will be an integral representation for 
the Dirac propagator in the non extreme Kerr-Newman metric and asymptotically at 
infinity, where the Kerr-Newman metric goes over in the Minkowski metric in oblate 
spheroidal coordinates. The construction of the propagators relies on Chandrasekhar 
separation of variables, whose completeness plays a crucial role in order to develop a 
reliable scattering theory. Their integral representations together with estimates for the 
asymptotic behavior of the solutions of the radial problem, arising in Chandrasekhar 
ansatz, will allow us to give explicit analytical expressions for wave operators classical 
at the event horizon and Dollard-modihed at infinity. This will be, to our knowledge, 
the first analytical result, since all previous works in this direction (see BED, which 
are based on the Mourre theory, treat only the problem of the existence of the above 
mentioned wave operators and of the asymptotic completeness for the Dirac equation 
in the Kerr-Newman metric. Our global strategy is the following: 

• in the first paper we prove the completeness of Chandrasekhar ansatz for the free 
Dirac equation in the Minkowski metric in oblate spheroidal coordinates and we 
derive an integral representation for the Dirac propagator in this framework; 

• in the second paper we give an integral representation for the Dirac propagator in 
the Kerr-Newman geometry in terms of the solutions of the radial and angular ODEs 
arising in Chandrasekhar separation of variables and we describe the behaviors 
of the radial solutions close to the event horizon and asymptotically at infinity, 
respectively; 

• in the last paper we treat first classical wave operators asymptotically at infinity 
and we compute them explicitly by means of the integral representations of the 
Dirac propagator obtained in the previous two works. In doing so, we get an 
analytical expression of the time-dependent logarithmic phase shift, that we need 
in order to construct Dollard-modified wave operators asymptotically away from 
the black hole. After implementation of this phase shift in the free dynamics 
we evaluate the Dollard-modified wave operators and obtain for them an integral 
representation. Then, we compute the wave operators close to the event horizon 
and prove the asymptotic completeness of the wave operators at the event horizon 
and asymptotically at infinity. 

The present article is structured as follows: In Section 1 we motivate our interest to 
derive an integral representation for the Dirac propagator in the Minkowski metric in 
oblate spheroidal coordinates. In Section 2 we introduce the Dirac operator in the Kerr- 
Newman metric. In Section 3 we find a suitable transformation, which decomposes 
it into the sum of an operator containing only derivatives respect to the variables t, 
r and tp and of an operator involving only derivatives respect to t, $ and ip. This 
transformation together with Chandrasekhar ansatz allows us in Section 4 to construct 
explicitly the generalized total angular momentum J for the Dirac equation in the 
Kerr-Newman background and to show that J is a symmetry operator. In Section 5 
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we restrict our attention to the free Dirac operator in the Minkowski metric in oblate 
spheroidal coordinates. In this framework we bring the Dirac equation into the form of 
a Schrodinger equation and we construct explicitly a positive scalar product respect to 
which the Hamilton operator is hermitian. In Section 6 the main results are Theorem 
6.2, where we prove the completeness of Chandrasekhar ansatz and Theorem 6.3, which 
gives the integral representation for the Dirac propagator in the Minkowski metric in 
oblate spheroidal coordinates. 

2. Preliminaries 

In the Newman-Penrose formalism pj in two-component spinor notation the Dirac 
equation coupled to a general gravitational field and a 4-vector field V is given in 
Planck units h c (■' 1 by i] 

(V aa' + ^ e ^AA')P A + W*Q*a' = 0 (1) 

(V Ay y - ieV AA >)Q A + i^P* A , = 0, (2) 

where V 4y y is the symbol for covariant differentiation, e is the charge or coupling 
constant of the Dirac particle to the vector held V, p* is the particle mass m e divided 
by y/2 and P A and Q A are the two-component spinors representing the wave function. 
Here, the asterisk used as a superscript denotes complex conjugation. Notice that the 
factor 2“ 2 in the definition of the mass /i* is due to the fact, that the Pauli matrices 

as defined in the Newman-Penrose formalism differ from their usual definitions by the 

factor y/2. Moreover, the vector potential enters with opposite signs in CO and Hi in 
order that gauge invariance be preserved, since the spinors in the above equations are 
related by complex conjugation. 

According to |3J, we denote with £ A a basis for the spinor space and with £*, A a basis 
for the conjugate spinor space. To the spinor basis we can associate at each point of 
the space-time a null tetrad (1, n, m, m*) obeying the orthogonality relations 1 • n = 1, 
m ■ m* = — 1 and 1 • m = 1 ■ m* = n • m = n ■ m* = 0. Furthermore, the covariant 
derivative of a spinor C, A can be expressed in terms of its components along the spinor 
basis £ A as follows 0) 0 


^X^AA B = d a A b + Kaa^ 


where T b , are the spin coefficients and 

caa x 


«■■'== a ==»‘£ 

(3) 

V :=* = »* A dw , :=s , :=m ,9_ 

(4) 


are the directional derivatives along 1, n, m, and m*. Following [4] and letting 
P° =: Fj, P 1 =: F 2 , Q* 1 ' =: G u Q*°' =: -G 2 , 
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equations © and m can be brought into the form 



CfoT = 0 




with 


^ — im e 

0 

A +/x* — 7 * +ieVt,n b 

— ( 6 * +/3* —r* +ieVj,m 

o D = 

0 

—im e 

— (<5+7t*— a* ) 

D+e*-p*+ieV b l b 

D+e—p+ieVbl b 

S* +7 T—at+ieVbm* b 

—im e 

0 


^ <5+/3— T+ieVi,m* b 

A+fi— 7+ze\4n b 

0 

—im e 


( 5 ) 


\ 
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( 6 ) 


and T = (iq, F 2 , G 1; In what follows, we consider the Dirac equation in the Kerr- 

Newman metric, i.e. in the presence of a rotating charged black hole. In Boyer-Lindquist 
coordinates (f, r, d, ip) with r>0, 0 < d < 7r, < ip < 2 ti the Kerr-Newman metric is 

given by d] 

sin 2 1 ?. 


ds 2 = — (dt — a sin 2 ddip) 2 — E ( —— + dd 2 

2.J \ ZA 


-[adf — (r~ + a 2 )d(^] (7) 


with 


E := E(r, d) = r 2 + a 2 cos 2 d, A := A(r) = r 2 — 2 Mr + a 2 + Q 2 , 

where M, a and Q are the mass, the angular momentum per unit mass and the charge 
of the black hole, respectively. In the non-extreme case M 2 > a 2 + Q 2 the function A 
has two distinct zeros, namely, 

r 0 = M — \J M 2 — a 2 — Q 2 and rq = M + \J M 2 — a 2 — Q 2 , 

the first one corresponding to the Cauchy horizon and the second to the event horizon, 
while in the extreme case M 2 = a 2 + Q 2 Cauchy horizon and event horizon coincide, 
since A has a double root at r\ = M. Throughout our work we will consider the case 
M 2 > a 2 + Q 2 and restrict our attention to the region r > rq outside the event horizon. 
Hence, A will be always positive. Notice that by setting M — Q — 0 in m the Kerr- 
Newman metric goes over into the Minkowski metric in oblate spheroidal coordinates 
(OSC), namely 

ds 2 = dt 2 — -dr 2 — A ^■^■d'd 2 + sin 2 ddip 2 ^ , (8) 

with A = r 2 + a 2 . In fact, by means of the coordinate transformation 

x = Vr 2 + a 2 sin d cos p, y = Vr 2 + a 2 sin d sin ip, z = r cos'd (9) 

(JBJ can be reduced to the Minkowski metric in Cartesian coordinates. Moreover, in 
the OSC the surfaces r =const are confocal ellipsoids, while the surfaces d =const 
are represented by hyperbolids of one sheet, confocal to the ellipsoids. Since 
sgn(^) =sgn(cos , d) with d € [0,7r], the surface d =const is actually given by a half¬ 
hyperboloid, truncated at its waist and laying in the half-space * ^ 0 according as 
d ^ At this point it is interesting to observe that at r = 0 the ellipsoids degenerate 
to the disk x 2 + y 2 = a 2 sin 2 d, whose boundary {(r, d, p)\r = 0,d = (and hence 
E = 0) corresponds to the set of points at which the Kerr-Newman metric becomes 
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truly singular. Furthermore, the surfaces p =const look like bent planes, which are 
approximately vertical for large r, but flatten out and become horizontal at the edge of 
the disk. For further details we refer to © and |H]. 

In order to write out explicitly © for the Kerr-Newman metric in Boyer-Lindquist 
coordinates, we make use of the Kinnersley tetrad ©, namely 


l b = 


r + a n a 
A ’ ’ ’A 


n 


= (r 2 + a 2 , —A, 0, a) 

1 


m b = 


ia sin d, 0 , 1 , i esc d) . 


y/2(r + ia cos'd) 

Then, the non vanishing spin coefficients in © are ©33 


( 10 ) 

(11) 

( 12 ) 


1 


P = 


r — ia cos d ’ 
ia 


7T = — —p 2 sind, 

V2 


P = 


T = 


p* cot d 

ia sin d 

y/2 £ ’ 


- 

f ‘ 2£ 


7 = P + 


M 


2£ 


and 


P 


a = 7T — p* = —j= ( iap sin d + - cot d 


where p* denotes the complex conjugate of p. Notice that while the vanishing of e is 
due to the particular choice of the tetrad m, (ED- m, the other spin coefficients k, 
a, A and v are zero according to the Goldberg-Sachs theorem HU. By means of © and 
© the directional derivatives along the tetrad are computed to be 


D — D 0l 


A = - 

6 = - 
5* =- 


ADp 
2£ : 

p* £ o 
p£ o 

7f’ 


n _ d 1 
D ° - Wr + A 


d 


d 


(r 2 + « )wt + fly- 


dr A 


at 

a 


(r 2 + a 2 )-— + a 


dp _ 

a 


at a<^_ 


c ■= — 

0 ■ ad 

a 


■ i . q a a 

i 1 a sm d— + esc d 7 — 

at op 


ad 


a 


a 


An := ttt + i [a sin d— + esc d—— 


at 


a</? 


From 


we know that the vector potential for a rotating charged black hole is 


V b = — {-Qr, 0,0, aQr sin 2 d) 

Z j 


and therefore 


A 2S 


bra 6 = hm* 6 = 0. 


Notice that the vector potential V is stationary and axisymmetric. Regarding the 
elements of the matrix operator Od, we find that 


A + p* — 7 * + ieV b n b 


A 

\/2E 


- r — M .eQr 
Do H-x-1 - * 


A 


A 
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— (5* + (3* — t* + ieVb r m b ) = p ^£ 0 + ^ co ^ ^ + * a P* s* 11 ^ 

— (5 + 7r* — a* + ieVbm b ) = p* ^£ 0 + - cot 

= \/2 fi 


D + e* -p* + ieV b l b 
D + e — p + feVj ,/ 6 
5* + 7 r — a + ieVbirt* b 
5 + [3 — r + ieVbm* b 
A + p — 7 + ieVbn b 


n * -eQr 
v iW 
= V2( Do -p-^ 


= - p ( £ 0 + - c °t $ 


= — p* ^£ 0 + - cot $ — -iap sin id 

A / ~ r — M eQr 
D 0 + —-— + « ^ 


V U ' A ' " A 

According to El we may replace the Dirac equation Op^(t, r, d, <p) = 0 by a modified 
but equivalent equation 

r, i?, <p) = 0 with W = r5'- 1 Ox,5, (13) 

where ip = S' -1 '!' = (F l5 F 2; Gi, G^) 7 " and T and S are non singular 4x4 matrices, whose 
elements may depend on the variables r and d. 

3. The decomposition of the operator VV 

In this section we show that it is possible to find non singular matrices S and T such 
that the operator W decomposes into the sum of an operator containing only derivatives 
respect to the variables t, r and <p and of an operator involving only derivatives respect 
to t, d and <p. There are several choices of S and T, which accomplish this. Here, we look 
for matrices S and T such that they simplify later on the transformation of equation 
(usd into its Schrodinger form. 


Lemma 3.1 Let r be positive with r > r\. There exist a non singular 4x4 matrix 


1*1 1*1 


S = cdiag(—p, 2 2 A 2 , 2 2 A 2 ,—p*), det(S) = 


DA 


with 0 / c€ C and a non singular 4x4 matrix 

T = dzap(p* -1 , —p* -1 , —p -1 , p -1 ), det(T) = E 

such that 

w = yv (t , 

with 







im e r 

0 

Vav + 

0 

0 

—im e r 

0 

v/AD 

\[KV_ 

0 

—im e r 

0 

0 

Vav+ 

0 

im e r 


(14) 


( 15 ) 










(16) 


W, 




( —am e cos i9 0 0 jC + \ 

0 am e cost? —0 

0 £ + — am e cos'd 0 

—£_ 0 0 arri e cos r d 


where 


_ d 1 
V± : ~ dr T A 


d 


d 


(r 2 + a 2 )— + a— - ieQr 


r d 1 . / . d d 

L± : = — H— cot v =F i [a sin v— + esc v— 
dv 2 \ at dip 


dt dip 
d_ 
dt 


Proof. Let us define the matrix 

S := diag (h(r, $), A(r, i3),a(r, i3), j(r, i3)) 

such that det(S') 7 ^ 0 and h, A, a, 7 at least C 1 ((r 1 ,+oo) x [0, 7 r]). The equation 
OdS'i/j = 0 gives rise to the system 

— vm e hF\ + (A + p* — 7 * + ieVbn b )aGi — (5* + (3* — r* + ieVpnri > )^G? i = 0 
—im e AF 2 — (5 + 7 r* — a* + ieVbm b )aG\ + (D + e* — p* + ieVbl b )^G 2 = 0 


—im e aG\ + (D + e — p + ieVbl b )h,Fi + (5* + 7r — a + feI4m* b )AF 2 = 0 

—im e r yG 2 + (6 + (3 — t + ieVbm* b )hFi + (A + p — 7 + ieVbn b )AF 2 = 0. 

Starting with (El, we find that 

(A + fi* — 7 * + ieVb‘n b )aGi = — 


(17) 

(18) 

(19) 

( 20 ) 


da r — M 

aV++ fr + ~^ a 1 Gl 

-(5* + (3* - t* + ieV b m b )^G 2 = p ( r yC+ + wx + iap* sin d 7 ) G 2 


A 


dd 


( 21 ) 

( 22 ) 


with 


_d_ _l_ 
+ dr A 
d 


d 


d 


(: r 2 + a 2 ) — + a— - ieQr 


£+ : = 7777 + - cot i9 — i \ a sin d— + esc i3—— . 


di3 


dt dip 
d 


d 


dt 


dip 


In order to obtain W as in the statement of this Lemma, we impose in m that 
da r — M 

^ + _ a _<t = 0 ’ 

whose solution is given by a(r,d) = a($)A ~2 with a(i3) 7 ^ 0 for every i3 G [ 0 , 7 r], while 
in Ea we require that 

d'y 

— + iap* sin t) 7 = 0 , 
dv 

whose solution is 7 (r,d) = —7 (r)p* with 7 (r) 7 ^ 0 for every r G (ri,+oo). Hence, (ET1) 
and (Ell simplify to 


(A + /i — / y H“ icv\)Ti jcCri — —--I^_|_Gri 


-(<£ + (3* - t* + ieV b m b )^G2 = - 


7 (r) 


2E 

£+G 2 . 


(23) 

(24) 
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Regarding cm we have 


— (5 + 7T* — a* + ieV b m b )aGi = ( cr£_ + — j G 


Ai 


dd 


ID + <:* - // + ieV/bG, = V2 Ae>_- h-/>*-) g 2 , 


(25) 

(26) 


where 


£>_ : = 


_i. 

9r A 


d 


d 


(r + a)—+ a-- ieQr 

at dp 


- 9 1 a ■ ( . a 9 .. 0 

L_ : = —— H— cot i7 + j asmd— + esc d Tr¬ 
ad 2 \ dt dp / 

By imposing a (id) = C 3 with some constant 0 ^ C 3 6 C and observing that y(r, id) = 
—c^p* with 0 7 ^ C 4 € C sets the last two terms in the parenthesis on the r.li.s. of 
equal to zero, cm cm (Esi and EEC) become 

(A + p* — 7 * + ieV b n b )aG\ = -—— D + G\ 

-{5* + (3*-r* + ieV b m b )^G 2 = - ^ C + G 2 

— (§ + 7i* — 01 * + ieV b m b )aGi = 

A 2 

(D + e* — p* + ieV b l b )^G 2 = - c A V2p*V_G 2 . 

Concerning csD, a short computation gives 


(D + e — p + ieV b l b )hFi = \/2 (hV_ + —- ph ) 


(h* + 7 T — a + ieV b m* b )AF 2 = — p ( A£ + + 


9/z 

9r 

9A 

9d 


(27) 

(28) 


and by similar reasonings as we did for (ESI) and cm we obtain h(r, id) = —ph(id) and 
A(?’, d) = c 2 A(r) with 0 ^ c 2 6 C. Finally, regarding (12UI) . we get 

- ( dh \ 

(5 + (3 — t + ieV b m* b )hFi = — p* \ hC _ + — — zapsind/zj Fi (29) 

(A + p - 7 + ieV b n b )AF 2 = - ^ (aD+ + ^ + F 2 . (30) 


2 S 


9r 


A 


By proceeding similarly as we did for ESD and ESI), we find that A(r, id) = c 2 A 2 and 
h(r, id) = —Cip with 0 7 ^ q 6 C. Hence, (1271) . (12H1) . (1221) and (13U1) become 

(H + e — p + ieV b l b )hFi = — c\\[2pT>_F\ 

(h* + 7 T — a + ieV b m* b )AF 2 = — c 2 pA _ 2 C + F 2 
(5 + (3 — r + ieV b m* b )hFi = 

(A + p - 7 + ieV b n b )AF 2 = - ^^^—V + F 2 . 


2 S 
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Then, the operator S l OS is computed to be 


( 


S~ l OS = 


—im e 

0 


vgc3P*A£ 


2ci 

C3 P* 


V 4 


£_ 


\/2ci p A 2 

c 3 

\ 

\ C 4 


£>_ 


-^c + 

r- C3 1 

\/2c2pA3 ^ 

2c4 


C2 

—im e 

0 


C4P_ 
--^4- 

* i 

V2c4p* 

C2 

0 


\ 


x>_ 


/ 


By choosing a non singular 4x4 matrix T as follows 
T = diagO* -1 , -p* -1 , -p~\ p -1 ), 
it can be checked that 

( im e i —nm e cos i9 0 p^A'D . 


V2c 3 


W = TS~ OS = 


C4 


V^C; 


0 

L VAV- 


-im e r-\-am e cos $ 


C2 


Cl 

V2C4 

C2 




Sav- 


C3 

\ 

k C/i 




—£ + 
C3 


v'AZV 


C4 2c 4 

If we impose the following conditions 
c 3 c 2 

C4 = 7! = 7! = C1=:C ’ 

we finally have 

^ im e r—am e cos $ 0 

0 —im e r-\-am e cost? 

vAxu £+ 

y —v'Ahy 

Tliis completes the proof. 


— im e r—am e cos $ 0 

0 im e f+am e cos / 


w = 


vWd+ 

— C- 

—im e r—am e cos $ 

0 


£+ \ 
Vav- 
0 

*m e r+am e cos $ J 


□ 


Clearly, the result obtained in Lemma fd.il is also valid in the extreme case, where A is 
simply (r — M) 2 . Without loss of generality we set c = 1. The next result, which holds 
for the Dirac equation in OSC, follows directly from Lemma Id. ll bv setting M = Q = 0. 


Corollary 3.2 Let r > 0. There exist a non singular 4x4 matrix 

S {osc ) = dia^(-p, 25 A"3,23 A"5,-p*), det(S {osc) ) = i 

Si A 

and a non singular 4x4 matrix 


such that 


with 


W 


(OSC) 

(t,r,<p) 


T = diag(p*~ 1 , —p* -1 , —p _1 , p _1 ), det(T) = S 


1/y(05C) _ yy(OSC) . yy(05C) 


/ 

im e r 

0 

\^v + 

0 

\ 


0 

—im e r 

0 

v/Azl 



\/m5_ 

0 

—im e r 

0 


V 

0 

VAv+ 

0 

im e r 

/ 


w, 


(OSC) 




(31) 


(32) 
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where is given by 1771) . A = r 2 + a 2 and 

~ d f d a d \ 

± ’ dr ^ \_ dt A dip) 

Notice that the angular part of the operator W( t)1 ? i¥ ,) remains unchanged, when we treat 
the Dirac equation in the Minkowski space-time with oblate spheroidal coordinates. 
Turning back to the Kerr-Newman metric, it can be easily checked that 

= 0, (33) 

which also implies 

[nww] = [nw,w] =o. 

An analogous commutator holds for the Dirac equation in OSC. Going back to 
Lemma o we observe that the matrix T splits into the sum 

r = T( r ) + T(^), (34) 

where T( r ) = diag(—r, r, r, —r) and T ^ = zadiag(— cos id, cos id, — cos id, cos id) satisfy the 
commutation relations 

[r (r) ,r w ] =o, (35) 

[r (r.), W(i 5l s» 5¥? )] = [r w ,w (t , r)V ,)] =o. (36) 


4. Construction of the symmetry operator J 


In this section we will make use of Chandrasekhar separation ansatz in order to construct 
a new operator J for the Dirac equation in the Kerr-Newman metric. We show that J 
commutes with the Dirac operator Od given by ©, thus being a symmetry operator for 
Od- Moreover, the separation constant turns out to be eigenvalue of J, whose physical 
interpretation is that one of the square root of the squared total angular momentum for 
a Dirac particle in the presence of a charged rotating black hole. 

Taking into account that the geometry of our problem is axial symmetric and considering 
stationary waves with energy u, the (t, ip ) dependence of the spinors A entering in (THfl) 
is given by 

ip(t, r, -d, cp) = e“ J V( fc+ 2) v3 A( r ’> $) (37) 


where k G Z is the azimuthal quantum number of the particle and A ( r >$) £ C 4 . 
Inserting (EZD in (Hi, it is straightforward to verify that ^(r, id) satisfies the following 
equation 


(w (r) + w w ) A = o, 

where 


/ 

im e r 

0 

Vav + 

0 

\ 


0 

—im e r 

0 




\/AT_ 

0 

—im e r 

0 


V 

0 

Vav + 

0 

im e r 

) 


(38) 


( 39 ) 
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and 


with 


W w = 


/ —am e cos'd 0 

0 am e cos d 

0 £+ 

0 


0 £+ \ 

0 

-am e cos d 0 

0 am e cos d J 


K(r ) = cn(r 2 + a 2 ) — eQr + ( A, + - ) a 


(40) 


- d .K (r) 

d \ ( 1 ' 

t± : = — +-cotd± Q(d), Q(d) = acjsind + f k + - ) escd. 


By defining -0(r, d) = (/i(r, d), / 2 (r, d), gi(r, d), g 2 (r, d)) T , equation (13%1) gives rise to the 
following systems of first order linear partial differential equations 


\/~Kv + gi + im e rfi + (£+<72 — ctm e cosd/i) = 0 
\/~Sx>_g 2 ~ *m e r/ 2 — (£_gi — am e cos d/ 2 ) = 0 
V'AP-Zi - im e rgi + (£+/ 2 - am e cosd^i) = 0 
V'AP+Zs + im e rg 2 - (£_/i - am e cos d# 2 ) = 0. 


Let us now define 


/i(r,d) = 7 i(r)di(d), / 2 (r, d) = cr 2 (r)r 2 (d), (41) 

di(r,d) = ai(r)/?i(d), d 2 (r, d) = e 2 (r)/x 2 (d). (42) 


If we substitute (52) and (B2D into the above system, it can be easily seen that the 
requirement of separability implies 

/3i(d) = di(d), e 2 (r) = 71 (r), /z 2 (d) = r 2 (d), aq (r) = cr 2 (r). 

According to the notation used in (2j, we set 

71 (r) = R-(r), 8 1 (d) = 5_ (d), a 2 (r) = R+ (r), r 2 (d) = 5+ (d). 

Hence, we get 

fi( r , d) = i 2 _(r)S_(d), / 2 (r,d) = i? + (r)S + (d), 

di (r, 0) = R+(r)S- (d), g 2 (r, d) = (r) S+ (d) 


and the system of partial differential equations become 

(\/A V + R + + im e rR-)S- + ( C + S + — am e cos dS'_)i?_ = 0 (43) 

(VKV-R- — zm e ri? + )S' + — (C-S- — am e cos 0S + )R + = 0 (44) 

(VA^fi. — im e rR + )S _ + (£ + S+ — am e cos dS-)R + = 0 (45) 

(\/AX) + i? + + im e rR-)S + — (C-S- — am e cos0S + )R- = 0. (46) 
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Introducing four separation constants Ai, • • •, A 4 as follows 

C + S + — am e cos dSC = — X±S_ (47) 

£_SC — am e cos $S + = + A 2 5' + (48) 

C + S + — am e cos dSC = — A . 3 S_ (49) 

C-S- — am e cos diS + = + A 4 S+, (50) 

we obtain from the set of equations (1751) . (ITTI) . (1771) and (1751) 

Vav + r + + im, e rR_ = Ai (51) 

\/A V + R + + im e rR_ = A 4 (52) 

y/A.f)_R- — im e rR + = A 2 R + (53) 

VAV_R_ — im. e rR + = A 3 i?+. (54) 

Clearly, the systems of equations la,-,® and (EH),-",® will be consistent if 
Ai = A 2 = A3 = A4 =: A. 

Notice that the decoupled equations give rise to the following two systems of linear first 
order differential equations, namely 


y/AV_ —im e r — A 
im e r — A \fKV + 


—C- A + am e cos 9 
A — cim e cos 6 C + 



= 0 . 


(55) 


(56) 


Starting from the systems of equations and m- ■ -,m we are now able 

to construct a new operator J. Moreover, we will show that such operator commutes 
with the Dirac operator Op. 

By setting Ai = A 2 = A 3 = A 4 =: A in (1771) .- • -. (15111) and (1571) .- • -. (1571) . multiplying 
(S3) by —rR_ and (ED by ia cos dSC, respectively, and summing together the resulting 
equations, we find that 

— iap* cos dy/~Kt> + gi + rp*C + g 2 = A/i (57) 

with p* = — (r + fa cost?) -1 . By means of the same method we obtain from the couples 
of equations 1 (1751) . (1521) 1 and (' (1551(1571) ) 


— iap* cosd\/A V_g 2 — rp* t^gi 

= a / 2 

(58) 

+iap cos'dVAV_f 1 + rpC+f 2 

= A. 9 , 

(59) 

+iap cos '^^/~Kv + f 2 — rpt^fi 

= A g 2 

(60) 


with p = —{r—ia cos d) 1 . Equations ( 1571 ) .- • •. ( 1601 ) give the entries of the matrix operator 
we are looking for, namely 

/ 0 0 — iap* cos 'dy/~K'D + rp*C + \ 

0 0 — rp*C _ —iap* cos'dy/AV- 

iap cos fiy/AV- rp£ + 0 0 

— rpt _ iap cos dy/AT> + 0 0 


J = 


\ 


/ 
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The operator j can be written in a more compact form as follows 

j = r _ 1 (r w w (r) - r (r) w w ). (ei) 

Clearly, a similar expression holds for the Dirac equation in OSC. 

Lemma 4.1 J = SJS _1 is a symmetry matrix operator for the Dirac operator Od in 
the Kerr-Newman metric, i.e. 

[O d , J] = 0. 

Proof. From m we have Od = ST 1 >VS 1 , which acts on the spinors T and we 
consider the operator J = SJS -1 with J given by m- ^ holds 

o d j = sr-'wr-'tTwWv - r ^w^s- 1 
jo d = sr-\r w w (r) - r (r) w w )r- 1 >v5- 1 

and therefore we obtain 

[O d , J] = ST^PS- 1 

with 

p := - r (r) w w ) - (r w w (r) - r (r) w w )r~ 1 >v. 

Making use of the decomposition W = W( r ) + W(#) and of the commutation relations 
we obtain 

p = w w (r (r) r -1 + - ^(r-'r^ + r^r- 1 )^^ 

= W(^)W( r ) - W(r)W(p) = 0 

where in the last line we employed <1531i. 

□ 

Concerning the physical meaning of J , it is interesting to observe that, when a = 0 the 

angular eigenfunctions S± (d) can be expressed in terms of the spin-weighted spherical 

harmonics Yf\ (see for instance IEMEI). According to m we will call J the square 
“^2 

root of the total squared angular momentum operator. In m it was proved that the 
separation constant A is the eigenvalue of the square root of the total squared angular 
momentum for the Dirac equation in any type-D vacuum space-time. Moreover, the 
occurrence of this operator and therefore the separability of the Dirac equation in the 
Kerr-Newman metric by means of Chandrasekhar ansatz arises from the presence of a 
Killing spinor field on the space-time under consideration. For a detailed description of 
the concept of Killing spinors we refer to Cl- 

5. Schrodinger form of the Dirac equation in OSC 


In order to bring the matrix equation 

yv {osc) tji = o 


( 62 ) 
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with W ( ° 5C) given by m into the form of a Schrodinger equation 
id t i> = H d ^, 


(63) 


we apply a method similar to that one used in HE]. Without risk of confusion we can 
simplify our notation by omitting to write explicitly the superscript ( OSC ) attached to 
the operator W. Starting with El, we bring the time derivatives on the l.h.s. of the 
equation and find that 


with 


- iTdti> = {W(r,v) + 


T = 


( 0 0 iV A 

0 0 —a sill'd 

— i\/~K —a sind 0 

^ —a siii-d z\/A 0 


—asind \ 
— 7a/A 
0 

0 / 


(64) 


and with W( ri¥ ,), formally given as in Corollary 15.21 where the operators V± and 

C± have been now replaced by 

~ dad 

± dr^Adtp 

~ d 1 . d 

= + x cot ^ =F «esc d—. 

ad 2 d<p 

Since det(T) = E 2 ^ 0 for r > 0 and d E [0,7r], T is non singular and we can multiply 
(El on both sides by — T 1 to obtain 

H D = ^S-(W {r , ip) +VV { o,v ) ), (65) 

where 02 is the Pauli matrix and 


asind ( 0 ^ 

’ ^ \ 0 -a 2 ) 


yV(r,<p) 

m e r ( 0 

1 2 

0 

) 

+ diag(- 

-£-, 





/ 

0 

0 

i 

0 \ 


( 

0 


am e cos d 


0 

0 

0 

i 



-M 

VA 


—i 

0 

0 

0 

+ 


0 



V 

0 

—7 

0 

0 J 


V 

0 


with 

£± = iD± and M.± = —=C± 

\'A 


such that 


£-) 

-M+ 

0 

0 

0 


0 

0 

0 

M- 


0 \ 

0 

M+ 

0 / 


E± = — £± and M.*± = 

Notice that the formal operator Hp acts on the spinors 1 /’ on the hypersurfaces t =const. 
Moreover, the matrix S entering in El is hermitian. For simplicity in the notation 
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we will omit the hat of the wave functions. It can be checked that for every 
</> G C£°([0,oo) x S ' 2 ) 4 the Hamiltonian H D is hermitian, i.e. formally self-adjoint with 
respect to the positive scalar product 

/»oo nl /» 27 r 


{ip \(j>) — / dr d( cosd) / d<p {ip\C<p) 

Jo J -1 Jo 

with inner product given by 

where the ip denotes the complex conjugated transposed spinor and 

a siii-d ( <72 0 

y/A \ 0 -cr 2 


C — H 4 + 


( 66 ) 

(67) 

( 68 ) 


In order to check the positivity of the scalar product defined by (ESj) and m it suffices 
to show that the matrix C entering in the inner product is positive definite. Indeed, the 
eigenvalues of C are given by 


Ai — A 2 — 1 


|a| sin d 

“7^ 


and A 3 = A 4 = 1 + 


|a| sind 

“ 7 ^ 


and since the following inequality holds for all r > 0 


lal sind 


< 


< 1 , 


yfK ~ yfK 

we can conclude that A, > 0 for every i — 1, • • • ,4. In what follows we consider the 
Hilbert space H = {C 4 , (■(•)} consisting of wave functions : [0, 00) x S 2 —> C 4 together 
with the scalar product 


6. Completeness of Chandrasekhar ansatz 

We begin with some preliminary observations. Energy, generalized squared angular 
momentum and the z-component of the total angular momentum form a set of 
commuting observables {Hp, J 2 , J z }. Moreover, the angular system (1561) can be brought 
in the so-called Dirac form [TO], namely 


US : = 


0 



sini!) 


am e cos d 
1 

au> sin d 


au sin d 


k+l 

sintf 

ain P cos d 


S = XS 


with S(d) = (S-(d), S + ( , d)) T and d G ( 0 , 7 r). According to pH] , in L 2 (( 0 , 7 t )) 2 
the angular operator U defined on D{U) = C^°(( 0 , 7 r )) 2 is essentially self-adjoint, 
its spectrum is discrete, non degenerate (i.e. simple) and depends smoothly on u. 
Therefore, its eigenvalues can be written as A j(u>) with j G Z and A j < A J+1 for 
every j G Z. Moreover, the functions S± (d) satisfy a generalized Heun equation 
PH] and become the well-known spin-weighted spherical harmonics by setting a = 0 . 
Furthermore, the functions e*( fc+5 ) l ’ c ’ are eigenfunctions of the z-component of the total 
angular momentum operator J z with eigenvalues — (k + |) with k G Z. Hence, we 
will label the states in the Hilbert space 7i by (J k j. In what follows, we consider the 









17 


free Dirac operator Hp given by m i n — D 2 ([0,oo) x S' 2 ) 4 . From Theorem 1.1 
fIT\ the operator H D defined on Cq°([0, oo) x S' 2 ) 4 is essentially self-adjoint and has a 
unique self-adjoint extension on the Sobolev space LF 1,2 ([0, oo) x S' 2 ) 4 . In addition, the 
spectrum of H Dl which we will denote by od, is purely absolutely continuous and given 
by &d = (—00 ,m e \ U [ra e ,+oo). Notice that, since Hp possesses a unique self-adjoint 
extension, when we will derive an integral representation for the Dirac propagator, there 
will be no need to impose Dirichlet boundary conditions or other boundary conditions 
on the spinors, as it was done in PI- In preparation of the proof for the completeness 
of Chandrasekhar ansatz we show in the next Lemma that 




silw 

st‘ + w 




(69) 


with k, j Eld is a complete orthonormal basis of L 2 (S 2 ) 2 . 


Lemma 6.1 For every k, j e Z the set </?)} with Yff 3 fd,ip) given by k (if A) is a 

complete orthonormal basis for L 2 (S 2 ) 2 . 


Proof. We begin by proving the orthonormality, i.e. we show that for every k', j', k, 
j £ Z it holds (Y$\Y* 3 ) S 2 = S kk >8^> , where is the usual scalar product on S 2 . 

Making use of (ESI), a direct computation gives 


(Y kj \Y k ' j ') s *=S kk ' [ dd smdS kj {V)S k ' j \tf), 

Jo 

where for simplicity in notation we omitted to write the subscript u attached to the 
angular eigenfunctions. In order to investigate the above integral we need to go back to 
the angular system (15 6|) . It is easy to see that S ^ J , S k f^ and ( S+) , (iSjM satisfy 
the following two systems of first order linear ODEs, namely 

+ £ + S'+ ,J + (Ay — am e cos 9)S k L =0, (70) 

—+ (Xj'+ am e cos 6) S’f-’ =0 (71) 


and 

+c + (s+y + (\j — am e cos 9)(s k _ j y = 0, (72) 

—+ (A j + am e cos 6) = 0, (73) 

where (•)* denotes complex conjugation and 

£± = Q(ti)± = ^coti?± Q(-d), Q(i?) = auosmd + fk' + ^ csci?. 


After multiplication of ((TUI) by sin d , of (ITU) by sin ti (jS+'j , of (ITU}) by sin r dS kj 

and of (ITU}) by sin , we consider the equations (ITUll sind — (ITUl) sin i)S k j = 0 

and (HU) sin $ \S+) — (ESI) sin 'dS'f 3 = 0, build their sum and integrate over $ from 0 
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to 7 r. Hence, we obtain 


2(A,/ - Aj) / d-d sintfS 1 ^)^ J (0) = 


l k 'j' - 


s k A c!_s k _ j -si j c- 


c\si j +s k _ j c + (s k J 


/ <M sin id 

Jo 

Since the angular eigenfunctions S±($) even vanish at $ = 0 and at id = it (see [2D]), it 
can be checked that the following relations hold 

+ J dd smdS k J C+ (s+ j Y = - J dd sind jtsA’ 

-J dd sindS^’ C_(s k jy = + J o d'J (s k jy jC + SA . 

By means of the above results we obtain that 

/*7T 

2 (Ay - A j) / did sin dS kj (d)S k ' j \'d) = 


S'?) s kj + 


U-CA S[ 


: k j 


did sin id 


and this completes the first part of the proof. Regarding the proof of the completeness, 
the method we use employs the Projection Theorem (Theorem II.3 l22]), i.e. we show 
that in L 2 (S 2 ) 2 the only element orthogonal to our orthonormal basis is the zero vector. 
Without loss of generality let us consider 


ip = 


( \ Ak+\)v 


V2n \ ^2 (d) 


with ipi e C£°((0,7r)) for z = 1,2. 


Then we have 

Ap\Y kj ) s 2 = J dx ((pi (x)S + (x) + ip\(x) 5_ (x)), 

where we made use of the transformation x = cos id. The following estimate holds 

\(Jp\Y k3 )sA < J dx (K(x)| |S+(x)| + W* 2 (x)\ |S1(x)|) 

< \Wi\\l 2 J ^dx \S + {x )\ 2 + \M 2 l 2 J ^dx |S_(x)| 2 

< d J dx (|5 + (x)| 2 + |R_(x)| 2 ) = d, 

where in the second line we used Holder inequality, in the third line the orthonormality 
condition for S±(d) and we defined d := max{||y>i||| 2 ,11y 2 11 2 2 }• Clearly, | (Lp\Y k ' J )s 2 | < 0 
if and only if d — 0. Since d = 0 implies ifi — ip 2 — 0 and the scalar product {-\-)s 2 is 
positive, the proof is completed. 

□ 


In addition from Theorem II.7 in El it follows that the Hilbert space L 2 (S 2 ) 2 is 
separable. We state now the Theorem on the completeness of Chandrasekhar ansatz. 
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Theorem 6.2 For every G C“([0, +oo) x S' 2 ) 4 it holds 
V>(0 ,r,d,p) = 


R k J,+(r)Y^(i), V ) 

\ RiuoY^ + («,<p) / 

with on = (— oo,m e ] U [m e ,+oo) and scalar product (•)•) defined by MIM . 


’ CJ D 


dai E Af) 

kJZzTj 


R k i = 


A 


Proof. Let us first recall that Lemma 6.1 implies that the space Cff(S 2 ) 2 is spanned 
by an orthonormal system {Yffifd, (^)}. We show now that it is possible to construct 
isometric operators 

W kJ : C 0 °°([0, +oo )) 2 — C 0 °°([0, +oo) x S 2 ) 4 , 

such that 

/ Rlh WWWv) \ 
flS+w 

R k j,+a)YpMv) 

\ RtU r ) Y F(A>p) / 

with A a positive dehnite hermitian matrix and 

Y*i(Av) = 

v Z7T 

Since the angular eigenfunctions Y k i are normalized, we have 

poo f‘ 

dr 4Si«£W) with pj = 


_ I Rtl M 
RiUa 


u Iln2([o,+oo)) 2 


Let us dehne a matrix 


is 2 


/ R kj Y k i \ 

Tokj v kj 
ll u ;,+ J lj,+ 
jokj ykj 
- rt CJ,+ 1 u- 

\ R%,Y% + ) 


. 2 A 
A 2 := — 


I 4 + 


a sind / cr 2 0 


VA 


0 


-u 2 


Clearly, it holds A 2 C = I 4 with C given by and A will be a positive dehnite 
hermitian matrix. Therefore, taking into account that the matrices C and A 2 commute, 
we obtain that 




f 27T' 


L 2 ([0,+oo)) 2 


dr I d( cos-d) / dip A 2 f) k J) 




'-1 


= A 2 ip k J) = = (W hJ (R k ‘)\W tJ (R kl )) 

= llW(«' y )lL( I 0 . + oo)x S =).- 

By means of the isometric operators Wk,j we can now introduce for every u G On an 
auxiliary separable Hilbert space tj(u;) as follows 

f)M = (J) f )kj with fj fcJ = B 4 j(C“([ 0 ,+oo)) 2 ). 

/cj'EZ 












20 


Moreover, the expansion theorem (Th.3.7 [23]) implies that 

t = E Wiwvy 

k,j£ Z 

Notice that f)(u;) is by no means a subspace of L 2 ([0,+oo) x S' 2 ) 4 with respect to 
the spatial measure, since the solutions R^±(r) of the radial system (ESD oscillate 
asymptotically for r —> +oo. To circumvent this problem we proceed as follows. First 
we identify the absolutely continuous part of with the operator H D , since the 
spectrum <jd is purely absolutely continuous. Then, lemma 10 (Ch.l §3.4 [23]) implies 
that H d is unitary equivalent to the operator of multiplication by u in L 2 ((7 d; da;) 4 . By 
definition (see Ch.l §5.1 [23j ) the Hilbert space 7 L C is decomposed into the direct integral 
/*© 

H= (74) 

J CT£) 

if there is given a unitary mapping T of the space hi onto S) but since the operator 
Hd is self-adjoint, the spectral representation theorem (Theorem 7.18, |24| i implies the 
existence of the isomorphism T and this completes the proof. 

□ 


The vector valued function ip u = (Fip)(u) is called the representative element of the 
element ip G hi in the decomposition fZ3D- The scalar product in S) will be introduced 
according to [23], as follows 

Jcr D 

where (•[•) denotes the scalar product given by (1331) . Furthermore, in S) we can introduce 
a norm by means of || • ||^ = (J doj || ■ ||jpa. 

In the next result we give the integral representation of the free Dirac propagator in 
Minkowski space in oblate spheroidal coordinates. 


Theorem 6.3 For every ip G C“([0,+oo) x S' 2 ) 4 and x = (r, $,</?) it holds 


ip{t,x) 


= e itHD ip(0,x) 



/cj'EZ 


with 


where on 



{r)Yi; + {d,<p) 

R k JAr)YF(A<p) 


\ RtU'FU^YP) / 

(—oo,m e ] U [m e , +oo) and the scalar product (•(•) is defined by J331) . 


(75) 


Proof. Since H D is self-adjoint in hi theorem 7.37 [23] implies that {U(t) = e ltHo \t G M} 
is a strongly continuous one-parameter unitary group with U{t)ip G C^([0, +oo) x S' 2 ) 4 
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for every i 6 I and every e C^([0,+oo) x S' 2 ) 4 . Making use of the version of the 
spectral theorem as given in 23| (Ch.l § 4-5), we have 

{(f>\e ltHD ip) — f du e lwt , (76) 

J C7D 

where {<j>u\ipu>) is given in terms of the resolvent of Ho as follows 

(</c|VO = 4~{<t>\E{v)i>) = ^-lim((j)\[R(uj - ie) - R{u + ie)]i/j) (77) 

CLUJ Am e—>0 

with E(uj ) the spectral family associated to Hr,. In addition Theorem 1.7 (Cli.10, §1 
[231 ) implies that {(f)J\^jJ) is absolutely continuous in awhile the unicity of the spectral 
family E(u) and the existence of a spectral family F(u>) = EE{uj)E~ l on S) follow 
directly from theorem 7.15 (23)- Making use of the relation (3) in (Ch.l, §4.2 ;23|) 
adapted to our case 

/ +oo 1 

—- dE(uj), 

■oo U- Z 


we get 


— -\R(u — ie) — R(uj + ie)] = — 
27 ri 7T 


r»+Od 


; dE(£u) 


(78) 


- (w - (7 ) 2 + e 2 

Since the integrand in (IZHD is bounded and integrable, when we take the limit for e —» 0 
in the above expression, we may apply Lebesgue dominated convergence theorem in 
order to take this limit under the integral sign. Taking into account that 
1 6 


lim ■ 

7T (uj — UJY + 6^ 


= 8{uj — ijj), 


we hnd that 


-—: liin[/?,(cj — ie) — R{uj + ie)) = Id. 
2m e—>o 


(79) 


Finally, (USD together with m and (1771) gives the desired result. 


□ 


In order to construct square integrable wave packets in the coordinate space, we go 
over in the Minkowski metric in Cartesian coordinates by means of the inverse of the 
coordinate transformations © and apply Theorem 1.8 in EH- Hence, we can choose 
H = L 2 ([0,+oo) x S' 2 ) 4 and normalize the wave packets according to ('0|'0)w = 1- 
The next Lemma describes locally uniformly in uj the asymptotics of the solutions of 
the radial system 


y/AV- —im e r — A \ / -R- N \_q 
im e r — A \J7s.V + J l R + J 


( 80 ) 


for r —>• +oo, where 


V± = — =F i 
dr 


uj + 


k±ii 

A 


A = r 2 + a 2 . 
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Lemma 6.4 Every non trivial solution R of \8tA) for |u>| > m e behaves asymptotically 
for r —> + oo like 


R(r) = 
with constants /“ 7 ^ 0 and 


R~(r) 

R+(r) 


cosh 0 sinh 0 
sinh 0 cosh 0 


e -inr[foo + 0{r~ l )\ 
e +iKr [ff° + 0{r~ l )\ 


1 f uj + m e 

Q = - log 

4 y uj — m e 


fuj_+rrie\ k _ e ( w )y/w 2 - m 2 , 
\lo - m e J 


where e(u) is a sign function such that e(uj) = +1 for uj > m e and e(uj) = — 1 for 
to < —m e . 

Proof. Let us first rewrite (18011 as follows 


R'(r) = V(r)R(r), V(r) = 
where * denotes complex conjugation and 


—zf 2 (r) (p{r) 
(p*(r) iQ(r) 


„ , . (k + 7 ) a A + im e r 

tt(r) = uj+ y A 2 , <p(r) = 


A 


C, C 2 
C 2 c i 


VA ' 

, R(r) = 


By means of the ansatz 

R(r) = AR(r), A = 
with constants C\ and C 2 such that Cf — Cf 7 ^ 0 , (IHT1) becomes 
R'(r) = T(r)R(r ), T(r) = i"V(r)i = 


A_(r) 

A4.fr) 


-Tu(r) -Ti 2 (r) 
Ti 2 (r) T n (r) 


where 


rn(r) = 


T 12 (r) = 


®(C? + Cf)fi(r) + CiCaf^fr) - <p(r)) 
Cl - C| 

2iC' 1 C' 2 fi(r) + Cf<p*(r) - C 2 V(r) 


Cl - Cl 

We £x now Cj and C 2 by requiring that 

lim Tn(r) — in, C\ — C\ = 1. 

r—>+00 

Taking into account that asymptotically for r —> + 00 , it results that 


T 11 M = 

condition (ESI) becomes 


i(Cl + C 2 2 )u - 2im e CiC 2 

Cl-Cl 


+ 0 -; , 


(C 2 + Cl)u — 2m e CiC 2 — k = 0. 

A simple calculation gives 

Ci = cosh 0, C 2 — sinh 0, 0 = - log 


uj + m e 
uj — m e 


(81) 


(82) 


(83) 
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With the help of (IH31) T n and T 12 can be rewritten as follows 


Tn(r) = - ( uQ(r) - m 2 e 

h\j 


Va 


T 12 (r) = 


A 

7 ^ 


,m e 


K 


+ i — fl(r) — a; 


Va )' 


Let us now introduce the decomposition T{r ) = A(r ) + B{r) with 
A(r ) = 


|a;2 m ev7 


—tt ~ im e - 
VA e k 


-4= + im e - 


1 - 


hA 


1 -Vk) 


^ 2 m eVA) 


and 


B(r) = 


i (k + |) 


-m P 


hi 


A 




Since A(r) and B{r ) are continuously differentiable on [0,+oo), 


p+oo 


dr\A (r)| < oo, 


r»+CX) 


dr\B(r)\ < oo 


and the eigenvalues Ai and A 2 of the matrix Aq dehned by 


A 0 := lim Air) = 

r —>+oo 


-IK 0 
0 IK 

are simple, then Theorem 11 in J2SJ, which is actually a simplified version of a theorem 
due to Levinson EH, implies that the equation R'(r) = T(r)R(r) has solutions such 
that for r —> +oo 

' e -iKr[foo + o{r~ l )\ 

e +iKr [f “ + 0(?’ -1 )] 

where /“ and /“ are the non zero components of the eigenvectors and £+ belonging 
to the eigenvalues — ik and +ik of the matrix A 0 , respectively. 

□ 


R(r) = 
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